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Abstract. Model atmospheres of isolated neutron stars 
with low magnetic field arc calculated with Compton scat- 
tering taking into account. Models with effective temper- 
atures 1, 3 and 5 MK, with two values of surface gravity 
(log g = 13.9 and 14.3), and different chemical compo- 
sitions are calculated. Radiation spectra computed with 
Compton scattering are softer than the computed with 
Thomson scattering at high energies (E> 5 keV) for hot 
(T ff > 1 MK) atmospheres with hydrogen-helium compo- 
sition. Compton scattering is more significant to hydrogen 
models with low surface gravity. The emergent spectra of 
the hottest (T g > 3 MK) model atmospheres can be de- 
scribed by diluted blackbody spectra with hardness fac- 
tors ~ 1.6 - 1.9. Compton scattering is less important for 
models with solar abundance of heavy elements. 



1. Introduction 

Relatively young neutron stars (NSs) with age < 10 6 yr 
are sufficiently hot (T c s ~ 1 MK) and can be observed 
as soft X-ray sources. Indeed, at present time the ther- 
mal radiation of few tens different kind of isolated NSs, 
from anomalous X-ray pulsars to millisecond pulsars are 
detected. The thermal spectra of these objects can be de- 
scribed by blackbody spectra with (color) temperatures 
from 40 to 700 eV (see, for example, Mereghetti et. al 
2002). 

The plasma envelope of a NS (if it exists) can be con- 
sidered as a NS atmosphere, and structure and emergent 
spectrum of this atmosphere can be calculated by using 
stellar model atmosphere methods (Mihalas 1978). Such 
modelling was performed by many scientific groups, be- 
ginning with Romani (1987), for NS model atmospheres 
without magnetic field as well as for models with strong 
(B > 10 12 G) magnetic field (see review by Zavlin & 
Pavlov 2002). These model spectra were used to fit the 
observed isolated NSs X-ray spectra (see review by Pavlov 
et al. 2002). 
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One of the important results of these works is as fol- 
lows. Emergent model spectra of light elements (hydrogen 
and helium) NS atmospheres with low magnetic field arc 
significantly harder than corresponding blackbody spec- 
tra. These elements are fully ionized in atmospheres with 
T er r > 1 MK. Therefore, the true opacity in these atmo- 
spheres (mainly due to free- free transitions) decreases with 
photon energy as E~ 3 . At high energies electron scattering 
is larger than true opacity and photons emitted deep in the 
atmosphere (where T > T c g) escape after few scatterings 
on electrons. In previous works, concerning isolated NS 
model atmospheres, coherent (Thomson) electron scatter- 
ing is considered. As a result, emergent spectra are very 
hard. But such a situation is very favorable to change the 
photon energy due to Compton down-scattering. 

It is well known that the Compton down-scattering de- 
termines the shape of emergent model spectra of hotter NS 
atmospheres with T e g ~ 20 MK and close to Eddington 
limit (London et al. 1986; Lapidus et al. 1986; Ebisuzaki 
1987). These model spectra describe the observed X-ray 
spectra of X-ray bursting NSs in Low-Mass X-ray Bina- 
ries (LMXBs), and they are close to diluted blackbody 
spectra with a hardness factor / c ~ 1.5 - 1.9 (London et 
al. 1986; Lapidus et al. 1986, Madej 1991; Pavlov et al. 
1991). But these model atmospheres with Compton scat- 
tering taken into account arc not calculated for relatively 
cool atmospheres with T c g < 10 MK. Therefore, the effect 
of Compton scattering on emergent spectra of isolated NS 
model atmospheres with T e s < 5 MK is not well known 
up to now. 

Here we present model atmospheres of isolated NSs 
with Compton scattering taken into consideration and in- 
vestigate the Compton effect on the emergent spectra of 
these atmospheres. 



2. Importance of Compton scattering 

First of all we consider the Compton scattering effect on 
emergent model spectra of isolated NS atmospheres qual- 
itatively. It is well known that in the non-relativistic ap- 
proximation (hu,kT e << m e c 2 ) the relative photon en- 
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Fig. 1. Comptonization parameter Zc omp vs. photon en- 
ergy for neutron star model atmospheres with different 
chemical composition. 

ergy lost due to a scattering event on a cool electron is: 



Each scattering event changes the relative photon energy 
by this value. It is clear that the Compton scattering effect 
can be significant, if the final photon energy change is 
comparable with the initial photon energy. Therefore, we 
can define the Comptonisation parameter Zc omp (see also 
Suleimanov ct al. 2006): 

Zcompt = max((T c *) 2 , r*), (2) 

m c c 

where max((T*) 2 , r*) is the number of scattering events 
the photon undergoes before escaping, r* is the Thomson 
optical depth, corresponding to the depth where escap- 
ing photons of a given frequency are emitted. We can ex- 
pect that Compton effects on emergent spectra of isolated 
NS model atmospheres are significant if the Comptoni- 
sation parameter approaches unity (Rybicki & Lightman 
1979). Because of this we compute Zc m P at different pho- 
ton energies (see Fig. 1) for hot NS model atmospheres 
with different chemical compositions. These models were 
computed by using the method described in the next sec- 
tion, with the Thomson electron scattering. It is seen from 
Fig. 1 that the Comptonisation parameter is larger (0.1 - 
1) at high photon energies (E > 4 — 5 keV) for H and 
He model atmospheres. Therefore, we can expect a signif- 
icant effect of Compton scattering on the emergent spec- 
tra of these models. On the other hand, Zcomp is low for 
the model with solar chemical composition of heavy ele- 
ments. The Compton scattering effect has to be weak on 
the emergent spectrum of this model. 

This qualitative analysis shows that Compton scatter- 
ing can be significant for light element model atmospheres 
of isolated NS we have investigated in more detail. 



3. Method of calculations 

We computed model atmospheres of isolated NSs assum- 
ing planar geometry using standard methods (Mihalas 
1978). 

The model atmosphere structure for a isolated NS 
with effective temperature T c g, surface gravity g, and 
given chemical composition, is described by the hydro- 
static equilibrium equation, the radiation transfer equa- 
tion, and the energy balance equation. These equations 
have to be completed by the equation of state, and also 
by the particle and charge conservation equations. We as- 
sume local thcrmodynamical equilibrium (LTE) in our cal- 
culations, so the number densities of all ionisation and 
excitation states of all elements have been calculated us- 
ing Boltzmann and Saha equations. We take into account 
pressure ionisation effects on the atomic populations using 
the occupation probability formalism (Hummer & Mihalas 
1988) as described by Hubeny et al. (1994). 

Compton scattering is taken into account in the ra- 
diation transfer equation using the Kompanccts operator 
(Kompaneets 1957; Zavlin & Shibanov 1991; Grebenev & 
Sunyaev 2002) . The energy balance equation also accounts 
for Compton scattering. 

For computing the model atmospheres we used a ver- 
sion of the computer code ATLAS (Kurucz 1970; Kurucz 
1993), modified to deal with high temperatures; see Ibrag- 
imov et al. (2003) for further details. This code was also 
modified to account for Compton scattering (Suleimanov 
& Poutanen 2006; Suleimanov et al. 2006). 

Our method of calculation was tested by compar- 
ing models for X-ray bursting neutron star atmospheres 
(Pavlov et al. 1991; Madej et al. 2004), with very good 
agreement. 



4. Results 

Using this method the set of hydrogen and helium NS 
model atmospheres with effective temperatures 1, 3, and 
5 MK and surface gravities log g = 13.9 and 14.3 were 
calculated. Models with Compton scattering and Thom- 
son scattering were computed for comparison. Part of the 
obtained results are presented in Figs. 2-4. 

The Compton effect is significant for spectra of hot 
{T e g > 3 MK) hydrogen model atmospheres at high ener- 
gies (Fig. 2). The hard emergent photons lost energy and 
heat the upper layers of the atmosphere due to interac- 
tions with electrons. As a result the high energy tails of 
the emergent spectra become close to Wien spectra, and 
chromosphere-like structures with temperatures close to 
color temperatures of the Wien spectra in the upper layers 
of the model atmospheres appear. Moreover, the overall 
emergent model spectra of high temperature atmospheres 
in first approximation can be presented as diluted black- 
body spectra with color temperatures that are close to 
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Fig. 2. Top panel: Emergent (unredshifted) spectra of 
pure H low gravity NS model atmospheres. Solid curves 
- with Compton effect, dashed curves - without Compton 
effect , dotted curves - blackbody spectra, thin dash-dotted 
curves - diluted blackbody spectra with hardness factors 
3.1, 1.87 and 1.65, for models with T cff = 1,3 and 5 MK. 
Bottom panel: Temperature structures of the correspond- 
ing model atmospheres. Effective and color temperatures 
are shown by dotted and dash-dotted lines. 



Wien tail color temperatures: 



7T 



B E (T C ), 



T c = f c T t 



Cff) 



(3) 



where / c is hardness factor. These results arc similar to 
those obtained for model atmospheres and emergent spec- 
tra of X-ray bursting NS in LMXBs (see Pavlov et al. 
1991; Madej et al. 2004). 

The Compton scattering effect on the emergent model 
spectra of high gravity atmospheres is less significant 
(Fig. 3). The reason is a relatively small contribution of 
electron scattering to the total opacity in high gravity at- 
mospheres compared to low gravity ones. 

The Compton scattering effect on helium model at- 
mospheres is also less significant than on hydrogen model 
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Fig. 3. Top panel: Emergent (unredshifted) spectra of 
H NS model atmospheres with different surface gravities 
(solid curves - high gravity models, dashed curves - low 
gravity models). For comparison the model spectra with- 
out Compton effect are shown for hottest high gravity 
model (dotted curve). Bottom panel: Temperature struc- 
tures of the corresponding model atmospheres. 

atmospheres with the same T e g and log g (Fig. 4) . The 
reason is the same as in the case of high gravity models. 
The contribution of electron scattering to the total opacity 
is less in the helium models. 

We also computed one isolated NS model atmosphere 
with solar metal abundances and T e g = 3 MK and 
log 5=14.3 (see Fig. 5). The model was calculated with 
Thomson and Compton scattering and we found that the 
Compton effect on the emergent spectrum is very small. 

5. Conclusions 

Emergent model spectra of hydrogen and helium NS at- 
mospheres with T ff > 1 MK are changed by the Comp- 
ton effect at high energies (E > 5 keV), and spectra of 
the hottest (T e g > 3 MK) model atmospheres can be de- 
scribed by diluted blackbody spectra with hardness fac- 
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Fig. 4. Top panel: Emergent (unredshifted) model spec- 
tra of pure He low gravity NS model atmospheres. For 
comparison the model spectra of pure H atmospheres arc 
shown by dashed curves. The model spectra of hottest 
atmospheres without Compton effect are shown by dash- 
dotted and dash-dot-dotted curves. Bottom panel: Tem- 
perature structures of the corresponding model atmo- 
spheres. 

tors ~ 1.6 - 1.9. But at the same time the spectral energy 
distribution of these models are not significantly changed 
at the maximum of the SED (at energies 1-3 keV), and 
effects on the color temperatures are not large. 

The Compton effect is most significant for hydrogen 
model atmospheres and to low gravity models. Emergent 
model spectra of NS atmospheres with solar metal abun- 
dances are changed by Compton effects very slightly. 
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Fig. 5. Emergent (unredshifted) model spectra of high 
gravity NS atmospheres with solar abundance of 15 most 
abundant heavy elements with (dashed curve) and with- 
out (solid curve) Compton scattering. The dotted curve is 
the corresponding blackbody spectrum. 
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